Quasiparticle properties of graphene in the presence of disorder 
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Abstract 

We calculate the quasiparticle properties of chiral two-dimensional Dirac electrons in graphene 
within the Landau Fermi Liquid scheme based on GW approximation in the presence of disorder. 
Disorder effects due to charged impurity scattering plays a crucial role in density dependence 
of quasiparticle quantities. Mode-coupling approach to scattering rate and self-energy in GW 
approximation for quasiparticle renormalized Fermi velocity and spin-antisymmetric Landau Fermi 
parameter incorporating the many-body interactions and the disorder effects show reduction of 
these quantities by 5 — 15 percent at available experimental charge carrier density region. 
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1. INTRODUCTION 



Graphene, a single atomic layer of graphite, has been successfully produced in experi- 
ment [l[ , which has resulted in intensive investigations on graphene-based structures because 
of fundamental physics interests and promising applications 2J. There are significant efforts 
to grow graphene epitaxially 3] by thermal decomposition of Silicon Carbide (SiC), or by 
vapor deposition of hydrocarbons on catalytic metallic surfaces which could later be etched 
away leaving graphene on an insulating substrate. 

The renormalized Fermi liquid parameters for both undoped and doped graphene have 
been firstly calculated by Das Sarma et al |4| focusing on the role of intraband density func- 
tions and on analytic weak-interaction expansions and they found that for doped graphene 
renormalization factor, and effective Fermi velocity show no deviation from the usual Fermi 
liquid behavior, and the Fermi liquid description is robust. On the other hand, with precisely 
zero doping, intrinsic graphene exhibits a quasiparticle lifetime linear in the excitation energy 
and a zero renormalization factor, indicating that the Fermi liquid description is marginal at 
the Dirac point. Next, the inelastic carrier lifetime in graphene {5] and many-body electron- 
electron quasiparticle properties in graphene have been correctly studied. 6] 

The quasiparticle energies and band gaps of graphene nanoribbons has been calculated 
using a first-principles Greens function approach within the GW approximation The 
electron self-energy, spectral function, and the band velocity renormalization due to phonon- 
mediated electron-electron interaction, has been calculated [8] and found that phonon- 
mediated electron-electron coupling has a large effect on the graphene band structure renor- 
malization and a reduction of the band velocity by 10 — 20 percent at the experimental 
doping level. 



aver of carbon are the effects 
9[ , the plasmon behavior and 



Other unusual many-body features of the single-atom-thick 
of electron-electron interactions on the ground-state properties 
essentially angle-resolved photoemission spectroscopy (ARPES) [lOj that follow from chiral 
band states and energy-momentum linear dispersion relation. 

The low energy quasiparticle excitations in graphene are linearly dispersing, described by 
Dirac cones at the edges of the first Brillouin zone. The Linear energy-momentum dispersion 
have been confirmed by recent observations ll|]. The slope of this linearity corresponds the 
Fermi velocity of chiral Dirac electrons in graphene and it plays essential role in the Landau 
Fermi liquid theory and has a direct connection to experimental measurement. ARPES 
is a powerful probe of spectral function, in two-dimensional crystals because it achieves 
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momentum resolution [12]. Recent experiments Q, Q, Q have reported ARPES spectra 
for single- layer graphene samples prepared by graphitizing the surface of SiC [la ]. Zhou 
et a/[l5| studied of the electronic structure of single layer epitaxial graphene with the main 
focus on the 7r bands by ARPES. The dispersion of the 7r bands shows overall agreement with 
conical dispersions of Dirac fermions. However, clear deviations from the conical dispersion 
were observed near the Dirac point. They examined the two possible explanations for the 
experimental data, many-body interactions and the opening of a gap at the Dirac point due 
to breaking of sublattice symmetry. Griineis et al\u\ have performed ARPES of graphite 
and compared the measured quasiparticle dispersion to ab initio calculations. They have 
found that the band dispersions are better described by the GW approximation, however 
around the K point that agreement becomes poorer and many-body approximation schemes 
beyond GW may improve result. From these experiments view, the renormalized velocity 
can directly be measured from the quasiparticle dispersion relation at long wavelength region 
near to the Dirac point. 

All the mentioned theoretical efforts calculating the quasiparticle quantities have been 
carried out for clean graphene. Since disorder is unavoidable in any material, there has 
been great interest in trying to understand how disorder affects the physics of electrons in 
material science specially here in graphene and its quasiparticle properties. 

Our aim in this work is to study the quasiparticle properties in the presence of electron- 
impurity and electron-electron interactions within the Landau Fermi Liquid theory where 
the concept of the model is applicable [18] to impurity systems and inhomogeneous electron 
liquids. It has been shown [l9( that the transport of undoped graphene cannot be described 
within the standard Boltzmann approach of metals, because the Dirac-like electronic ex- 
citations have infinite Compton wavelength violating the assumptions for the validity of 
Boltzmann transport. Accordingly the proper way to compute the transport and quasipar- 
ticle properties are treating the impurities in a self-consistent way. For this purpose, we use 
the self-consistent theory of Gotze [20|] to calculate the scattering rate, self-energy and the 
quasiparticle properties of the system at the level of random-phase approximation (RPA) 
including disorder effects. Implementing the impurity effects in our calculation, is in the 
same way of our earlier work on both chiral two-dimensional electron systems (C2DES) 21] 
where we have calculated the ground-state properties and thermodynamics of graphene and 
conventional two-dimensional electron gas (2DEG). 
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In this work, we consider the charged impurity and the surface-roughness potentials which 
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are always present in currently available graphene samples [23|, |24j . It has been demonstrated 
that a short-range scattering potential is irrelevant for electronic properties of graphene 25] 
and impurity scattering in graphene is dominated by charged impurities and it is therefore 
strongly dependent on the carrier charge density 26|. We have used the same method 
22|, |27| to investigate some properties of the conventional 2DEG. In this paper, we find that 
there are significant differences between quasiparticle properties of C2DES and conventional 
2DEG due to disorder effects. 

The rest of this paper is organized as follows. In Sec. II, we introduce the models for 
self-consistent calculation of impurities effect. We then outline the calculation of self-energy 
and quasiparticle quantities. Section III contains our numerical calculations. We conclude 
in Sec. IV with a brief summary. 



2. THEORETICAL MODEL 

We consider a system of C2DES interacting via the Coulomb potential e 2 /er and its 
Fourier transform v q = 2ne 2 /(eq) where e is the background dielectric constant at zero 
temperature in the absence of electron-phonon and spin-orbit interactions. The massless 
Dirac band Hamiltonian of graphene can be written as 7i = v r (o~i k\ + 02 k 2 ), where r = ±1 
for the inequivalent K and K' valleys at which n and it* bands touch, fcj is an envelope 
function momentum operator, v is the Fermi velocity and Oi is a Pauli matrix which acts on 
the sublattice pseudospin degree-of-freedom. 

Our results for C2DES are based on the GW-RPA in which the self-energy is expanded 
to first order in the dynamically screened Coulomb interaction W: 

s' m=—oo 
'1 + SS' COS (0 k ,k+q) 



X 



} Gi? ) (k + q,^ n + ^ m ), (1) 

where s = + for electron-doped systems and s = — for hole-doped systems, ft = l/(k B T), 
u n = (2n + 1)tt/ ft is a fermionic Matsubara frequency, G^j\k, iu) is noninteracting Green's 
function and the effective electron-electron interaction in GW approximation is given by 

W (q, ,0, F) = v q + v q Y^^—^^ , (2) 

where x^(q, zfi, /i, T) is the dynamical polarizability of noninteracting C2DES in which 
chemical potential and the strength of damping are represented by \i and T, respectively. In 
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the C2DES, dielectric-function contributions from intraband and interband excitations are 
subtly interrelated. The two contributions must be included on an equal footing in order to 
describe C2DES physics correctly. 

A central quantity in the theoretical formu 



ation of the many-body effects in Dirac 
fermions is the dynamical polarizability 9|, |28|, |29|] tensor yWf q, iQ, n ^ 0). This is de- 
fined through the one-body noninteracting Green's functions. |30j To describe the properties 
of C2DES we define a dimensionless coupling constant a gr = ge 2 /veh where g = g v g s = 4 is 
the valley and spin degeneracy. 

The effect of disorder is to dampen the charge-density fluctuations and results to modify 
the dynamical polarizability tensor. Within the relaxation time approximation the modified 
X^(q, /i, T) is given by 31] 



x (0) (q,^,/i,r) 



X (°)(q,^ + T» 



n+r 



X<°) (q,ifl+ir» 
X (°)(q) 



(3) 



We consider long-ranged charged impurity scattering and surface roughness as the main 
sources of disorder. The latter mechanism also known as ripples comes either from thermal 
fluctuations or interaction with the substrate. [32|] The disorder averaged surface roughness 
(ripples) potential is modeled as (\U sur f(q)\ 2 ) = 7rA 2 h 2 (27re 2 n/e) 2 e~ q2A2 ^ where h and A 
are parameters describing fluctuations in the height and width, respectively. We can use 
the experimental results of Meyer al. 23j who estimate A ~ 10 nm and h ~ 0.5 nm. It is 
important to point out that there are other models to take into account the surface-roughness 
potential. The effect of bending of the graphene sheet has been studied by Kim and Castro 
Neto 33j. This model has two main effects, firstly the decrease of the distance between 
carbon atoms and secondly a rotation of the p z orbitals. Due to bending the electrons 
are subject to a potential which depends on the structure of the graphene sheet. Another 
possible model is described by Katsnelson and Geim 25j considering the change of in-plane 
displacements and out-of-plane displacements due to the local curvature of a graphene sheet. 
Consequently, the change of the atomic displacements results to change in nearest-neighbour 
hopping parameters which is equivalent to the appearance of a random gauge field described 
by a vector potential. These different models need to be implemented in our scheme and to 
be checked numerically to assess their validity in comparison to the available measurements. 

The charged disorder potential is taken to be {\Ui mp (q)\ 2 ) = riiV 2 e~ 2qd in which is the 
density of impurities and d is the setback distance from the graphene sheet. 

We modify the mode-coupling approximation introduced by G6tze[20] for C2DES to 
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express the total scattering rate in terms of the screened disorder potentials 

~<| U imp (q) | 2 > 

q 

(\Usurf(q)\ 2 Y 



v F k F 

ll = - > 

2hnA ^ 



+ 



e"(q) 



£ 2 (q) 



(4) 



l + ir^ (q,ir)/x°(q) ' 
where A is the sample area, e(q) = 1 — v q x^(q) is the static screening function and 
the relaxation function for electrons scattering from disorder is given as <po(q., iT) = 

[x (0) (q^r,/i)-x (0) (q)] /»t . 

Since the scattering rate T depends on the relaxation function ipo(q, iT), which itself 
is determined by the disorder included response function, the above equation needs to be 
solved self-consistently to yield eventually the scattering rate as a function of the coupling 
constant. 

In Eq. (pP) the sum runs over all the bosonic Matsubara frequencies fl m = 2mn/(3. The 
factor in square brackets in Eq. §T§ , which depends on the angle #k,k+ q between k and k + q, 
captures the dependence of Coulomb scattering on the relative chirality ss' of the interacting 
electrons. The Green's function Gs (k, iuS) = l/[iuj — £ s (k)] describes the free propagation of 
states with wavevector k, Dirac energy £ s (k) = svk — /i (relative to the chemical potential) 
and chirality s — ±. The first and second terms in Eq. ([2j) are responsible respectively 
for the exchange interaction with the occupied Fermi sea (including the negative energy 
component), and for the interaction with particle-hole and collective virtual fluctuations. 
After continuation from imaginary to real frequencies, iuo — ► u + irj, the quasiparticle weight 
factor Z and the renormalized Fermi velocity can be expressed [18| in terms of the wavevector 
and frequency derivatives of the retarted self-energy S^(k,a;,r) evaluated at the Fermi 
surface (k = k F ) and at the quasiparticle pole u = £+(k): 

1 



- l-5^eEf(k,o;,r)| fc=fep>a)=0 ' 

and 

..* 1 4- dA- 1 ^ ^^ r , et (k,c;,r)L , 



(5) 



(6) 



v Z- 1 
On the other hand, we find that the quasiparticle renormalized Fermi velocity within 
the on-shell approximation is given by v*/v = 1 + (w) -1 9fc3?eE^ t (k, u>, T) | fc _ fcF w _ + 
9 w 3?eE^ t (k, to, T) \ k _ k ^_ - Evidently, the on-shell approximation is a valid expression for 
the quasiparticle renormalized Fermi velocity only in the weak coupling limit where a gr 
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FIG. 1: (Color online) The renormalized Fermi velocity scaled by Fermi velocity of a noninteracting 
clean system as a function of the coupling constant a gr for cut-off value A = k c jkp = 100 for 
different charged impurity densities value. In the inset, the renormalized Fermi velocity scaled by 
noninteracting one as a function of the coupling constant a gr for cut-off value A = 100 calculated 
within on-shell approximation. 



becomes small enough. The physical reasons are discussed in Ref. 34| in the context of 
conventional 2DEG. 

Following some standard manipulations [3] the self-energy can be expressed in a form 
convenient for numerical evaluation, as the sum of a contribution from the interaction of 
quasiparticles at the Fermi energy, the residue contribution S res , and a contribution from 
interactions with quasiparticles far from the Fermi energy and via both exchange and virtual 
fluctuations, the line contribution E lme . In the zero-temperature limit 

d 2 q v q \1 + s'cos (6> k ,k+q) 



k, w ,r 



(27r)2 £ (q,cu-Mk + q),r) 
x [9(o;-£ s ,(k+q))-e(-6,(k + q))] 



(7) 
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and 



d 2 q 

(2k) 



2 U 1 



1 + s' COS ($k,k+q) 



6'(k + q) 



, 2vre(q,«fi,r) [cu-^(k + q)] 2 + ft 2 

(8) 

where e(q, zO, T) = 1 — v q x^ \(l,iQ, fi,T) is the dynamic dielectric function in RPA level. 
Note that at the Fermi energy <9fc£+ s (k, u, T) vanishes, and c^E^^k, u, V) involves an inte- 
gral over interactions on the Fermi surface that are statically screened. These expressions 
differ from the corresponding 2DEG expressions because of the relative chirality dependence 
of the Coulomb matrix elements, because of the linear dispersion of the bare quasiparticle 
energies, and most importantly because of the fast short-wavelength density fluctuations 
produced by the interband contribution to x (l>^,r). we can introduce an ultraviolet 
cutoff for the wave vector integrals k c = Akp which is the order of the inverse lattice spac- 
ing and A is dimensionless quantity. Fermi momentum is related to density as given by 
k F = (AKn/g) l l 2 . 
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FIG. 2: (Color online) The quasiparticle weight factor Z as a function of the coupling constant a gr 
for cut-off value A = k c jkp = 100 for different charged impurity densities. 
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FIG. 3: (Color online) The spin- antisymmetric Landau Fermi parameter Fq as a function of the 
coupling constant a gr for cut-off value A = k c jkp = 100. In the insert, the spin-antisymmetric 
Landau Fermi parameter Ffi as a function of the coupling constant a gr calculated within the 
on-shell renormalized Fermi velocity approximation. 



In order to have in connection with experiments, we would like to calculate the spin- 
antisymmetric Landau Fermi parameter, Fq which is related to the spin-susceptibility, in 
which Xo/Xs — v*(l + Fq)/v. For this purpose, we obtain the Landau interaction function 
/k k' are defined by the following functional derivative [ls| 



./k.k' 



where the occupation numbers n+ ak in conduction band is denned as 0(kp CT — k). The 
spin-antisymmetric Landau Fermi interaction function f a ((p) is defined by — /kk')- 

After some straightforward algebraic calculations, f a {(p) on the Fermi surface is given by 



5n 



(o) 



(9) 



,cr',k' 



fM 



1 + COS (0k'-k,k' 



IF(k'-k,0) 



(10) 
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where <p being the angle between k and k' and |k| = |k'| = kp- Furthermore, Fq calculated 
using the following equation 

V*7l l n Z7T 



3. NUMERICAL RESULTS 

In this section we present our numerical calculations for quasiparticle properties of 
graphene in the presence of impurities which are modeled as mentioned above. The quasi- 
particle properties are calculated by using the theoretical models described above and the 
impurities dependence of quasiparticle quantities are calculated. In all numerical calcula- 
tions we consider d = 0.5 nm. Electron density is taken to be 1 x 10 12 cm~ 2 for Figs. 1-3. 

Increasing disorder (increasing rij or decreasing d for charge- disorder potential or increas- 
ing h for surface roughness potential) decrease the (<7> Q A*> F) as the scattering rate T 
gets bigger. Thus, decreasing f2, /x, T) (or increasing correlation effects) results in a 

stronger disorder potential. Despite r increases with increasing a gr , apparently it turns to a 



saturation limit and does not diverge 21(. This behavior is quite different than what is seen 
in conventional 2DEG [22] when the many-body effects influence the scattering rate through 
the local-field factor. In the conventional 2DEG system, at a critical level of disorder this 
nonlinear feedback causes T to increase rapidly and diverge, which is taken as an indication 
for the localization of carriers. However, in graphene, our calculations show that the T does 
not diverge therefore impurities cannot localize carriers and we have a weakly localized sys- 



tem in the presence of impurities compatible with experimental observations 35j . We have 



found through our calculations [2JJ that T increases with increasing rii/n as a function of 



Otg r . 



The Fermi liquid phenomenology of Dirac electrons in graphene [6, lOj and conventional 



0,Q 



2DEG 



have the same structure, since both systems are isotropic and have a single 
circular Fermi surface. The strength of interaction effects in a conventional 2DEG increases 
with decreasing carrier density. At low densities, the quasiparticle weight factor Z is small, 



the velocity is suppressed [34J ]. the charge compressibility changes sign from positive to 
negative 22], and the spin-susceptibility is strongly enhanced 34j. These effects emerge 
from an interplay between exchange interactions and quantum fluctuations of charge and 
spin in the 2DEG. 

In the C2DEG clean graphene, it was shown jf], 0, Q that interaction effects also become 
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noticeable with decreasing density, although more slowly, the quasiparticle weight Z tends 
to larger values, that the velocity is enhanced rather than suppressed, and that the influence 
of interactions on the compressibility and the spin-susceptibility changes sign. These quali- 
tative differences are due to exchange interactions between electrons near the Fermi surface 
and electrons in the negative energy sea and to interband contributions to Dirac electrons 
from charge and spin fluctuations. 
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FIG. 4: (Color online) The renormalized Fermi velocity scaled by Fermi velocity of a noninteracting 
clean system, if jv as a function of the electron charge density n in unite of 10 12 cm -2 at a gr = 1. 
In the inset, the scattering rate T as a function of electron charge density n at a gr = 1 for density 
impurities, 10 10 ,10 11 and 10 12 cm~ 2 , respectively. 

We have calculated the renormalized Fermi velocity scaled by Fermi velocity of noninter- 
acting clean system as a function of a gr in the presence of disorder as shown in Figure 1. 
It is found that the disorder effects become more appreciable at large coupling constants, 
within the mode coupling approximation. The renormalized Fermi velocity becomes smaller 
due to disorder effect about 10 — 15 percent at available experimental a gr value, namely 
smaller than 2 ( a gr ~ 1 or 2 is a typical value thought to apply to graphene sheets on the 
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surface of a (Si/Si0 2 or SiC) or SiC>2 substrate, respectively). Since the scattering rate T 
tends to a constant at large a gr , hence the effect of disorder for v*/v is fixed at large a gr . In 
the insert of figure 1, we have shown the results of renormalized Fermi velocity calculated 
within the on-shell approximation. As we pointed out the on-shell approximation is valid 
at small value of a gr . The behavior of v* shows some novel physics, which is qualitatively 
different from the physics known in the conventional 2DEG. In the conventional 2DEG, 
effective mass decreases by increasing impurity at small coupling constant, namely r s < 1 
and then enhances by increasing impurity effects at large r s value. 27 ] 

Figure 2 shows the quasiparticle weight factor Z as a function of a gr . It appears that 
the disorder effects become more appreciable at large coupling constant and decreases by 
increasing the disorder effects. Significantly, the dJR.eY7+ l : (k, u, T) does not change by impu- 
rity, however <9^9ffeE^ ne (k, u, T) has negative values which decreases by increasing impurity 
density. Moreover, the <9/ c !fteE^ Tie (k, u, T) has positive values which increases slightly by 
increasing impurity density. Accordingly, the quasiparticle weight factor Z decreases by 
increasing the impurities effect. Importantly, this behavior is in contrast with 2DEG where 
Z increases by increasing impurity effects. 

Figure 3 shows the spin-antisymmetric Landau Fermi parameter as a function of a gr by 
using the renormalized Fermi velocity calculated in Eq.(9). As it is clear from Eq.(14), 
this quantity is proportional to inverse renormalized Fermi density and its absolute value 
increases by increasing density. Consequently, the spin-susceptibility decreases by increasing 
impurity density as charge compressibility behaves 21J. In the inset, F£ as a function of 
coupling constant is shown by using the on-shell renormalized Fermi velocity. 

In Fig. 4 we show our theoretical predictions for the renormalized Fermi velocity of 
doped graphene scaled by Fermi velocity of noninteracting clean system as a function of 
density n in units of 10 12 cm -2 at a gr = 1. In connection with experiments, the gate 
voltage V, and electron density are related through the gate capacitance per area, C by 
V = en/C. Therefore, the results of v*/v as a function of n or V could be directly tested by 
experiments. For definiteness we take A = k c /kp to be such that 7r(A/cj?) 2 = (2tt) 2 /Ao, where 
A = 3V3ao/ 2 is the 

area of the unit cell in the honeycomb lattice, with ao — 1.42 A the 
carbon-carbon distance. With this choice A ~ (gn^ 1 x 10 2 , where n is the 

electron density in units of 10 12 cm -2 . We have shown the related scattering rate as a 
function of electron charge density in the inset of figure. 

As it is clear from Fig. 4 the renormalized Fermi velocity decreases by increasing impurity 
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density. By increasing the interaction effects, i.e., decreasing the density our theoretical 
effective Fermi velocity results decrease. Furthermore, we have found that at low electron 
density regime by increasing the impurity density, the renormalized Fermi velocity decreases 
sharply at = 10 10 cm~ 2 and reaches to a constant. Similar behavior occurs at larger 
impurity densities as well. Consequently, the impurity dependence of the renormalized 
Fermi density is quite nontrivial and nonmonotonic and this could be directly experimentally 
tested. Although, it is known that application of RPA in very low density in graphene is 
questionable, however it would be interesting to explore this effect by going beyond of RPA 
or using a density functional type model calculation by incorporating correctly the physics of 
linear dispersion relation, chirality and the exchange-correlation energy. It would be useful 
to carry out further experimental work at low electron density region to assess the role 
played by correlation and impurities effects. 



4. CONCLUSION 



We have studied the quasiparticle properties of the chiral electron graphene sheet within 
the Landau Fermi theory based on the random phase approximation incorporating the impu- 
rities in the system. Our approach is based on a self-consistent calculation between impurity 
effect and many-body electron-electron interaction. We have used a model surface roughness 
potential together with the charged disorder potential in the system. We find that there are 
significant differences between quasiparticle properties of C2DES and conventional 2DEG 
due to disorder effects. Our calculations of renormalized Fermi velocity show a reduction 
about 5 — 15 percent at experimental available density region. 

We remark that in a very small density region, the system is highly inhomogeneous and 
the effect of the impurities are very essential. A model going beyond the RPA is necessary to 
account for increasing correlation effects at low density. Our finding at low electron density 
would be verified by experiments. 
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